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ABSTRACT 

We use recent data obtained by three (OSSE, BATSE, and COMPTEL) of four instruments on board the Comp- 
ton Gamma Ray Observatory, to construct a model of Cyg X-l which describes its emission in a broad energy 
range from soft X-rays to MeV 7-rays self-consistently. The 7-ray emission is interpreted to be the result of 
Comptonization, bremsstrahlung, and positron annihilation in a hot optically thin and spatially extended region 
surrounding the whole accretion disk. For the X-ray emission a standard corona-disk model is applied. We show 
that the Cyg X-l spectrum accumulated by the CGRO instruments during a ^4 year time period between 1991 
and 1995, as well as the HEAO-3 71 and 72 spectra can be well represented by our model. The derived parameters 
match the observational results obtained from X-ray measurements. 

Subject headings: elementary particles — gamma rays: theory — plasmas — radiation mechanisms: non-thermal 
— scattering — stars: individual (Cyg X-l) 
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1. INTRODUCTION 

One of the brightest sources in the low-energy 7-ray sky, Cyg 
X- 1 , has been exten sively studied during the last three decades 
since its discovery (|Bowyer et al. 1965j, for a review see O da 
1977, Liang & Nolan 1984| ). It is a high-mass binary system 
(HDE 226868) with an orbital period of 5.6 days consisting 
of a blue supergiant an d presumably a black hole (BH) with a 
mass in excess of 5M Q ( Dolan 1992|). The separ ation of the two 



components is as 4 x 10 lz cm (IBeall et al. 1984). A periodicity 



of 294 d found in X-ray and optical light curves is thought to be 
related to pr ecession of the accretion disk ( Priedhorsky, Terrell, 
& Holt 1983, |Kempet al. 1983| ). 

The X-ray flux of Cyg X-l vari es on all observe d timescales 
down to a few milliseconds (e.g., Cui et al. 1997), but the av- 
erage flux exhibits roughly a two-modal behaviour. Most of its 
time Cyg X-l spends in a so-called 'low' state where the soft 
X-ray luminosity (2-10 keV) is low. The low-state spectrum is 
hard and can be described by a power-law with a photon index 
of ~ 1.7 in the 10-150 keV energy band. There are occasional 
periods of 'high' state emission, in which the spectrum consists 
of a relatively stable soft blackbody component and a weak and 
variable hard power-law component. Remarkable is the anticor- 
relation bet ween the soft and hard X-ray components (L iang & 
Nolan 1984), which is clearly seen during the transition phases 
between the two states. 

Cyg X-l is believed to be powered by accretion through an 
accretion disk. Its X-ray spectrum indicates the existence of a 
hot X-ray emitting and a cold reflecting gas. The soft black- 
body component is thought to consist of thermal emission from 
an o ptically thick and cool accretion disk (Shakura & Sunyaev 
1973, pringle 1981] |Batuciriska-Church et al. 1995|) . The hard 
X-ray part [J> 10 keV) with a break at ~ 150 keV has been 
attributed to thermal emission of the accreting matter Comp- 
tonized by a hot corona with temperature from tens to hun - 
dred keV ( |Sunyaev & Titarchuk 1980j |Liang & Nolan 1984 - 
A broad hump peaking at ~ 20 keV (Done et al. 1992), an 
iron Ka e mission line at ~ 6.2 keV with an equivalent w idth 



100 eV flBarr, White, & Page 1985| , |Kitamoto et al. 1990i see 



also |5bisawa et al. 199*6 and refe r ences therein ), and a strong 
iron K-edge (e.g., see|lnoue 1989| [Tanaka 1991 , Ebisawa et al. 
1992,1996) have been interpreted as signatures of Compton re- 
flection of hard X-rays off cold accreting material. 

In addition, there have also been sporadic reports of a hard 
spectral component extending into the MeV region. The most 
famous one was the so-called 'M eV bump' observ ed at a 5cr 
level during the HEAO-3 mission (Ling et al. 1987). For a dis- 



cussion of the pre-CGRO data and 7-ray emission mechanisms 
see, e.g., areview by Owens & McConnell (1992). The COMP- 
TEL spectrum accumulated over 15 weeks of real observation 
time during the 1991 -95 time period shows si gnificant emission 
out to several MeV (McConnell et al. 1997), which, however, 



remained always by more than an order of magnitude below the 
MeV bump reported from the HEAO-3 mission. 

The annihilation line search provided only tentative (1.9(t) 
evidence for a weak 511 ke V line with a flux of (4.4 ± 2.4) x 
10~ 4 photons cm' 2 s -1 (|Ling & Wheaton 1989[ ). Recent 
OSSE observations ( ^hlips et al. 1996| ) resulted only in upper 
limits with values of < 7 x 10~° cm s _1 for a narrow 511 
keV line and < 2 x 10~ 4 cm~ 2 s _1 for a broad feature at 51 1 
keV. 

Although an unified view for the X-ray spectra of BH candi- 
dates and their spectral states has yet to be constructed, the qual- 
itative picture seems to be quite clear. Current popular models 
include an o ptically thick disk component, a hot Comp tonizing 



region (e.g., Haardt et al. 1993 Gierlihski et al. 1997), and/or 



an advection-dominated accretion flow (e.g., Abramowicz et al 



1995, Harayan&TTTO 



and references therein). The spectral 
changes are probably governed by the mass accretion rate (e.g., 
|Chen et al. 1995| , foin, McClintock, & Narayan 1997| ). 

This picture, however, provides no explanation for the ob- 
served 7-ray emission (e.g., McConnell et al. 1997). The hard 
MeV tail can not be explained by standard Compton models be- 
cause they predict fluxes which are too small at MeV energies, 
and thus another mechanism is required. The models devel- 
oped so far connect the 7-ray emission with a compact hot core 
(~ 400 keV or more) in the innermost part of the accretion disk, 
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Table 1 
Luminosity of Cyg X- 1 . 



Energy band 


Luminosity, 10 36 erg/s 


> 0.02 MeV 


26 


0.02-0.2 MeV 


20.5 


0.2-1 MeV 


4.8 


> lMeV 


0.6 







which e mits via bremsstrahlung, Compton scattering, a nd anni- 
hilation ( |Liang & Dermer 1988 , [Skibo & Dermer 1995[ ), or with 
7T° productio n due to collisions of ions wit h nearly virial tem- 
perat ure (e.g., Kolykhalov & Sunyaev 1979], J ourdain & Roques 
1994). Li, Kusunose & Liang (1996) have shown that stochastic 
particle acceleration via wave-particle resonant interactions in 
plasmas (~ 100 keV) around the BH could provide a suprather- 
mal electron population, and is able to reproduce the hard state 
MeV tail. The possibility of Comptonization in the relativistic 
gas inflow near the BH horizon has been discussed by Titarchuk 
&Zannias(1998). 

We use the recent data obtained by three of four instruments 
aboard CGRO to construct a model of Cyg X-l, which de- 
scribes its em ission in a wide energy range from soft X-ra ys to 
MeV 7-rays ( [vloskalenko, Collmar, & Schonfelder I99l\ ). In- 
stead of a compact (pair-dominated) 7-ray emitting region, we 
consider an optically thin and spatially extended one surround- 
ing the whole accretion disk. It produces 7-rays via Comp- 
tonization, bremsstrahlung and positron annihilation. For the 
X-ray emission the corona-disk model is retained. 

In section 2 we discuss the combined OSSE-BATSE- 
COMPTEL spectrum of Cyg X- 1 . Our model and the inferred 
results are described in sections 3—4, and the implications are 
discussed in section 5. The applied formalism is given in the 
Appendix. 

2. OBSERVATIONS 

Since its launch in 1991 Cyg X-l has been observed by 
CGRO several times. The time averaged COMPTEL spectrum 
based on all observations between the CGRO Phases 1 and 3 
(April 19 91 to November 1994) is shown in Fig. |i] (M cConnell 
et al. 1997) together wi th the nearly con temporaneous spectrum 
derived from BATSE ( |Ling et al. 1997 ). The thick solid curve 
shows the best fit to the OSSE spectrum (0.06-1 MeV) for all 
obse rvations between April 1991 and May 1995 (P hlips et al. 
1996). The best-fit parameters for a power-law model with an 
exponential cutoff are a power-law photon index of T — 1 .39, a 
cutoff energy E c = 158 keV, a normalization intensity of 0.470 
photons cm" 2 s" 1 MeV" 1 at 0.1 MeV. The COMPTEL data 
provide evidence for a hard power-law tail extending up to at 
least 3 MeV. 

Although the OSSE and BATSE spectra have similar shapes, 
their intensity normalizations are different by a factor of ^2 
(Fig. [I]). The discrepancy is largest at the highest energies 
around 1 MeV. The COMPTEL measurements lie in between 
OSSE and BATSE. Although, there is no way of deducing the 
exact spectral shape in this region, the total spectrum is prob- 
ably smooth, without bumps, which is illustrated by the three 
individual spectra. Possible reasons for this discrepancy have 



been discussed by McConnell et al. (1997). For our further 
analysis we will use the combined BATSE-COMPTEL spec- 
trum. 

Table [I] lists the average luminosities of Cyg X- 1 for vari- 
ous energy bands. The values are derived from the combined 
BATSE-COMPTEL spectrum assuming a source distance of 
2.5 kpc. The total luminosity is between 1% and 10% of the 
Eddington luminosity, 



L Edd ee AttGM 



m p c 
ax 



1.25 x 10^ erg/s 



M 



IOMt 



(1) 



3. THE MODEL 



The existence of a compact pair-dominated core around the 
BH in Cyg X-l is unlikely in view of the CGRO observations. 
The sign ature of such a core would be a bump (L iang & Der- 
mer 1988, Liang 1990^ similar to the one reported by HEAO-3. 
However, no evidence for such a bump wa s detected by CGRO 
(|Phlips et al. 1994|McConnell et al. 1997|). In addition, the lu- 



minosity of Cyg X-l above ~ 0.5 MeV, though small, would 
substantially exceed the Eddington luminosity for pairs, which 
is ~ 2000 times lower than that for a hydrogen plasma. On the 
other hand, the hard MeV tail observed by COMPTEL can not 
be explained by Comptonization in a corona of kT ~ 100 keV 
and therefore another mechanism is required. 

Our study is an attempt to extend the 'standard' disk-corona 
model, which has been shown to work quite well at X -ray 
energies (e.g., Gierlihski et al. 1997, Dove et al. 1997), by 
including the processes of 7-ray emission. We investigat e 
the proton-dominated optically thin solution (Svensson 1984), 



G = kT/m e c <j 1, where the 7-ray emission is attributed 
to a spatially extended cloud surrounding the whole accretion 
disk (Fig. fy, the outer corona, which emits via bremsstrahlung, 
Comptonization, and positron annihilation. We concentrate on 
the hard X-ray to 7-ray part of the spectrum, and thus we in- 
clude into consideration the above-mentioned processes as well 
as Comptonization of the soft X-ray disk emission in the 'stan- 
dard' inner corona. The optical depth of the outer corona has to 
be small enough to avoid effective reprocessing of the emission 
from the disk and the inner corona. 

The soft X-rays consist of two components, the local black- 
body emission from the disk plus the reflected spectrum. At 
energies above ss30 keV the former is negligible and the later 
is only of minor importance. Therefore, we neglect both com- 
ponents at the moment and leave the detailed spectral modelling 
until the discrepancy in the intensity normalization of the OSSE 
and BATSE spectra has been resolved. However, the effective 
temperature of the soft excess is used in calculations of Comp- 
tonization in the inner and outer coronae, and the estimate of the 
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Fig. 1. — Diamonds show the average spectrum of Cyg X-l 
based on all COM PTEL CGRO Phase 1-3 observations (M c- 
Connell et al. 19977 an d crosses represe nt the almost contem- 
poraneous BATSE data ( Ling et al. 1997 ). The solid curve rep- 
resents the best fit to the time-averag ed OSSE spectrum con- 
taining 122 days of observation time. (Phlips et al. 1996). 



Fig. 2, 



A schematic view illustrating the model. 



total soft X-ray luminosity is provided to match the observed 
value. 

Our idea is that the electrons in the outer corona (which is 
optically very thin) are relativistic. The heating mechanism is 
not specified, but people usuall y refer to stochastic acceleration 
( Li, Kusunose & Liang 1996 ), MHD turbulence in the inner 
corona ( [Li & Miller 1997 ), and plasma instabilities in magne- 
tized advection- dominated accretion flows (B isnovatyi-Kogan 
& Lovelace 1997). These mechanisms are likely to heat mainly 
electrons and so can provide a population of energetic parti- 
cles. We further show (see Discussion) that the mean free path 
of these electrons in the outer corona is of the same order as 
its size. Because there is no mechanism to confine energetic 
electrons (and pairs if they exist), except for reasons of charge 
conservation, they can move freely inside the outer corona pro- 
viding the same temperature for the whole plasma volume. Ad- 
ditionally, the electron cooling in a thermal plasma at low num- 
ber density and small optical depth is not very efficient. 

We do not consider the process of ir° production in 
pp-collisions. Although it could be important at a few 
Schwarzschild radii (where the energy of protons is nearly 
virial), it is unimportant at tens to hundreds of Schwarzschild 
radii which is the characteristic size of the outer corona. The 
protons in the accreting flow far from the BH horizon should 
be cold, since the gravitational forces are quite weak there and 
thus the viscous heating in the disk is negligible. The energy 
transfer due to the Coulomb coupling with the hot electrons is 
also not efficient. 

3.1. The Fitting Parameters 

A set of eight fitting parameters was chosen: kTi, Ti, and 
kT a , t , the temperature and optical depth of the inner (/) and 
the outer (o) coronae, which are assumed to be spheres, L* oft , 
the luminosity of the disk which is effectively Comptonized by 
the inner corona, L so f t , the total effective soft X-ray luminosity 
of the central source illuminating the outer corona, R, the outer 



corona radius, and, Z = n + /n p , the positron-to-proton ratio in 
it. 

The formulae to calculate the bremsstrahlung, annihililation, 
and Comptonization emissivities are given in the Appendix. 
The accretion disk spectrum, which is further reprocessed by 
the inner and outer coronae, was taken to be monoenergetic 
with an energy Eq = l-SkT^ corresponding to the maximum 
of the Planck distribution, where kT^ — 0.13 keV is the ef- 
fecti ve temperature of the soft excess (B aluciriska-Church et al. 
1995). 

The bremsstrahlung and annihilation photon fluxes from the 
outer corona are proportional to R 3 ninj, where nij_ are the 



number densities of the plasma particles (see eqs. [A1],[A3]) 



Thus, if the annihilation contributes significantly, there is a con- 
tinuum of solutions given by the set of equations 



-/?7~ Z ( 1 ~\~ Z) 

R 3 n-n + = R 3 n 2 p Z{l + Z) = T ° 1 „ „J = . oust. 



2zy 



t = o~TRrip(l + 2Z) — const, 



(2) 



where kT is fixed, n p is the proton number density, gt is the 
Thomson cross section, and Z(l +Z)/(1 + 2Z) 2 varies slowly 
for Z > 0.5 (therefore, the fitting procedure is not very sensi- 
tive to this parameter). If only a negligible positron fraction is 
present, the continuum of solutions is defined by 



t = Rn p aT = const, 



(3) 



where kT Q is fixed, R < R max , and i? max is fixed from the 
fitting procedure. 

4. RESULTS 

The observed spectra of Cyg X-l are shown in Figs. [}] and^ 
together with our model calculations. The Comptonized spec- 
trum from the outer corona is only shown up to 3 MeV because 
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Fig. 3. — Upper panel shows the calculated Cyg X-l spec- 
trum together with the data points, which are the same as in 
Fig. D Central and l ower panels: the HEAO-3 72, and 71 spec- 
tra ( |Gng et al. 1987 ). In all panels the thin solid lines represent 
our model fit for the parameter sets I. The individual spectral 
components are the annihilation line (dotted line), ee-, e + e~- 
, ep-bremsstrahlung (dash-dot), and the Comptonized spectra 
from the inner and outer coronae (dashed lines). The Comp- 
tonized spectra from the outer corona are shown up to 3 MeV, 
up to which the approximation used agrees with Monte Carlo 
simulations, and where also significant data points are available. 
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Fig. 4. — The same as in Fig. ^, but for parameter sets I (see 
Table |). 
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Table 2 

The 'best-fit' model parameters. 





CGRO Phase 1-3 


HbAU-3: 


72 -state 


7i 


■state 


Pqi*q IT"] pfpr^ 

l al allic LCI a 


I 


n 


T 
1 


11 


r a 
i 


rr 
11 


Soft X-ray luminosity, L so f t (10 36 erg s _1 ) 




O.U 


10.6 


10.7 


9.8 


7.9 


/-corona temperature, fcT; (keV) 


76.7 


79.7 


95 


94.9 




93.0 


/-corona optical depth, Tj 


2.39 


2.23 


1.41 


1.42 




1.44 


L s * oft ,10 36 ergs" 1 


0.73 


0.84 


1.96 


1.95 




0.51 


o-corona temperature, kT a (keV) 


396 


436 


450 


448 


346 


361 


o-corona optical depth, r D 


0.06 


0.05 


0.056 


0.056 


0.12 


0.10 


o-corona radius, i? (10 8 cm) b 




^ioo 


^ioo 


150 


391 


812 


Positron-to-proton ratio, Z h 











1.0 


1.0 


0.5 


Proton number density, n p (10 10 cm _3 ) b 


k, 900 


>750 


>840 


187 


154 


93 


Accretion disk radius, Rj (10 s cm) 








1 


1 


1 


511 keV line flux, I a (1CT 5 photons cm -2 s^ 1 ) 











0.18 


0.15 


0.04 


xl 


4.0 


3.9 


1.4 


1.4 


0.9 


0.9 



a The inner corona is small or even absent at all 
b For R, n p , Z dependence see eqs. (Q), ([3]) 



= The narrow annihilation line flux from the disk (eq. [ A4]) as calculated for the given Rj 
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obser ved luminosity of the soft excess (B alucihska- 
1. 1995, for a distance of 2.5 kpc). This value agrees 
covering factor <^ 0.2 obtained by Doveetal. (1997) 
insistent Monte Carlo modelling of the corona-disk 
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e produce the observed broad-band X-ray spectrum 
(pierlinski et al. 1997|, [bove et al. 1997^. 
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licture is supported by X-ray observations. 

analysis of source temperature (defined from 
brems strahlung model) vs. 45-140 keV inten- 



et al. 1996) showed that the source temperature and 

130 



the intensify vary only within a limited range: 
keV, and 



170 



0.07 — 0.12 photons cm~ 2 s _i , with few low- 
temperaturfe - low-amplitude exceptions. A similar behaviour 
of the best- fit bremsstrahlung temperature vs. hard X-ray lumi- 
nosity (40- 200 keV) has been found by Kuznetsov et al. (1997) 
from the aialysis of the entire dataset of the Granat/SIGMA 
observatioi s of Cyg X-l collected between 1990 and 1994. 
The soft X-ray (< 10 keV) lu minosity of Cyg X-l is on the 
1.5 x 10 36 erg/s (e.g., |Liang & Nolan 1984], E bisawa 



Duri ng the HEAO-3 71, 72 states it was even lower 



1987). Taking into account that for the hard X-ray 



photons the Compt onization efficiency in t he hot plasma drops 
substantial y (e.g., fr-lua & Titarchuk 1995 ) and the number of 



2 The value of xt can not be used for the likelihood criterion estimates here mainly because of the uncertainty in the relative normalization of the OSSE, BATSE, 
and COMPTEL data. 
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Fig. 5. — The transport cross section for electrons in a hydro- 
gen plasma vs. the Lorentz factor of a particle. The individual 
lines correspond to the plasma temperatures (from top to bot- 
tom): 9 = 0.2, 0.6, 0.8, 1.0. For comparison the Thomson 
cross section gt is also shown. 



Fig. 6. — Upper panel: The total cooling rate due to the Comp- 
ton scattering as function of the plasma temperature for several 
values of the optical depth (eq. |B4[ ). Lower panel: The cool- 
ing rates due to the electron bremsstrahlung w ee + w ep (ER, 
NR) and Coulomb coupling with cold protons vs. plasma tem- 
perature (Z — 0, see Appendix). For the Compton scattering 
shown are the average cooling rates calculated for R = 10 10 
cm, where the line styles and the optical depths correspond to 
these in the upper panel. 
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photons decreases as well, the values of L so ft ~ 10 37 erg/s we 
obtain match the observational results. 

No pairs are required to reproduce the spectrum of Cyg X- 
1 in its normal state. If one takes an annihilation line flux of 
I a « 4.4 x 1CT 4 photons cm" 2 s" 1 (|Ling & Wheaton 1989j) 



t a {0 ~ 1) w 400 s 



10 12 cm" 3 



(4) 



in the 71 state, the accretion disk radius is estimated to be 
Rd ~ 1.7 x 10 9 cm (eq. [\4], set I). The upper limit allowed 
by optical measureme nts is Rd w 6 x 10 9 cm (M/10Mq) 
(Liang & Nolan 1984), while the effective radius of the soft X- 
ray emitting regio n of the disk is ~ 4.6 x 10 7 cm (B aluciriska- 
Church et al. 1995). 

Our calculations show that the presented model is consistent 
with the available observations of the Cyg X-l system, and is 
able to reproduce the observed spectra well. A more detailed 
study, however, would require a solution of the discrepancy in 
the intensity normalization between the OSSE and BATSE data 
and further Monte Carlo modelling. 



5. DISCUSSION 

We have calculated the radiation from Cyg X-l self- 
consistently assuming that the hot optically thin outer corona 
exists. A mechanism of its maintenance was not specified in 
our model (so it is not totally self-consistent), but the energy 
required to maintain such a corona is quite small and could be 
provided by a turbulent mechanism, stochastic particle accel- 
eration, and/or dif fusion of high energy electro ns from the in- 
ner disk (e.g., see fr^i, Kusunose & Liang 1996 and references 
therein). A relevant example is the solar corona of ~ 10 6 K 
(though its energetic contents is low) compared to 6000 K of 
the Sun's effective temperature; but a direct scaling to a BH is 
not appropriate here. In this section we discuss the physical 
conditions in the outer corona, i.e. diffusion of electrons and 
the cooling mechanisms, while do not touch its origin. 

To treate the diffusion of energetic electrons in the outer 
corona we consider the transport cross section for ee- 
scattering. This allows us (i) to exclude unimportant scattering 
at small angles dominating in Coulomb interactions, and (ii) 
provides us with a correct estimate of the typical cross section 
since the ep-collisions in a hot plasma are of minor importance 
compared to ee-collisions. 

The transport cross section a tr (see Appendix) for electrons 
in a hydrogen plasma of = 0.2, 0.6, 0.8, and 1.0 is shown in 
Fig. |[ Although, the value of ot r for particles with a Lorentz 
factor 7 ~ 2 is much larger than the Thomson cross section, the 
corresponding mean free path of electrons is close to the radius 
of the outer corona. This allows electrons to pass freely and 
therefore provide the same temperature for the whole plasma 
volume. Positrons, if produced somewhere, can also homoge- 
neously fi ll the plasma vo lume. The annihilation time scale is 
given by ( |Svensson 1982 ) 



1 



-(1 + 26 2 ln _1 [1.12e + 1.3]), 



For the parameters listed in Table || t a is of the order of hundred 
of seconds. 

The relevant cooling rates for electrons in a pure hydrogen 
plasma Z = (see Appendix) are shown in Fig. g (lower 
panel). The rates are divided by n 2 , the Coulomb logarithm 
was taken as a constant of In A = 20. For comparison we 
show the average Compton energy losses per unit volume (di- 
vided by n 2 ), wc m rip 2 = W Cm rip 2 x 3/(4ttR 3 ) = W Cm x 
3cr|/(47ri?r 2 ) (see eq. [|[), calculated for L soft = 10 37 erg/s, 
t = 0.04, 0.05, 0.06, and R = 10 10 cm. Clearly the average 
Comptonization losses (for r D fixed) depend on the radius of the 
outer corona and for 9 > 0.4 and the parameters adopted sub- 
stantially dominate bremsstrahlung losses and losses due to the 
Coulomb coupling with cold protons. On the other hand, their 
total value is not too high, of the order of I/ so ft (upper panel), 
which is about 10-20% of the total luminosity. 

The average value of the Compton energy losses of an elec- 
tron in the outer corona can be estimated as dE/dt = Wcm/N, 
where N = 4ttt R 2 /3ctt is the total number of electrons in the 
outer corona. Taking the corresponding numerical values (Ta- 
ble @), W Cm » L soft ~ 10 37 erg/s, r - 0.5, R = 10 10 cm, 
one can obtain dE/dt w 200 keV/s. The appropriate timescale 
for an electron of 7 ^ 2 is few seconds which is long compared 
to R/c <1 sec the particle needs to cross the outer corona. 



6. CONCLUSION 

The data obtained recently by the CGRO instruments allow 
us to construct a model of Cyg X-l which describes its emission 
from soft X-rays to MeV 7-rays self-consistently. This model is 
based on the suggestion that the 7-ray emitting region is a hot 
optically thin and spatially extended proton-dominated cloud, 
the outer corona. The emission mechanisms are bremsstrah- 
lung, Comptonization, and positron annihilation. For X-rays a 
standard corona-disk model is applied. 

The CGRO spectrum of Cyg X-l accumulated over a ~4 
years period between 1991 and 1995, as well as the HEAO-3 
71, and 72 spectra can be well represented by our model. The 
derived parameters match also the basic results of the X-ray 
observations. A fine tuning of the model would require further 
Monte Carlo simulations and more accurate spectral measure- 
ments. In this respect, the solution of the discrepancy between 
the OSSE and BATSE normalizations would be of particular 
importance. 

We thank the referee for useful comments. Discussions with 
R.Narayan, L.Titarchuk, and M.Gilfanov are greatly acknowl- 
edged. We are particularly grateful to M.McConnell for pro- 
viding us with the combined spectra of Cyg X-l prior to publi- 
cation, and E.Churazov for Monte Carlo simulations of Comp- 
tonization in 9 ~ 1, t k 0.1 — 0.05 plasma. 



APPENDIX 

A. RADIATION FROM A THERMAL PLASMA 

For a thermal plasma consisting of electrons, positrons and protons at mildly relativistic temperatures (kT < m e c 2 ), the main 
radiation processes are bremsstrahlung, electron-positron annihilation, and Compton scattering. 
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A. 1 . Bremsstrahlung 

The bremsstrahlung emissivities, the number of photons emitted per unit time, per unit volume, and per unit energy interval, can 
be represented by the form 

e - B /e 



Sij (e, kT) oc riiirij 



-G i3 {e,kT), 



(Al) 



where e = E/m e c 2 is the dimensionless photon energy, O = kT/m e c 2 is the dimensionless plasma temperature, and n^j with 
i = {e _ , e + }, j = {e~, e + ,p} are the corresponding number densities. Accurate numerical fits for the Gaunt factors G ee (e, kT) 
and G e + e - (e, kT) in an appropriate energy range have been given by Stepney & Gui lbert (1983) and Haug (19 87), respectively. The 
ep-bremsstrahlung emissivity can be calculated by the one-fold integration (e.g., see Stepney & Guilbert 1983). The approximations 
of the Gaunt factors Gy (e, kT) have also been constructed by Skibo et al. (1995). 



A. 2. Comptonization 

To calculate the effect of Compton scattering in a medium of kT ~ 100 keV we follow the model by Sunyaev & Titarchuk (1980) 
with corrections made by Titarchuk (1994) and Hua & Titarchuk (1995). The total number of photons emerging from the plasma 
cloud per unit energy interval and per unit time is given by 



F{E,kT) = 



F„(x,xo)L soft 



x «1, x <C x, 



(A2) 



where x = E/kT, xq = Eo/kT, E is the photon energy, Eq is the energy of soft photons injected into the plasma, L so ft is the 
lumi nosity of the soft photon source, and F v (x, xp) is the emergent spectrum represented by the Green function (H ua & Titarchuk 
1995). 

The results of the Hua & Titarchuk (1995) model are generally in a good agreem ent with Monte Carlo simulations except at high 
temperatures, O ~ 1, and small optical depth, r ~ 0.1 — 0.05 ( [Skibo et al. 1995 ). However, it still provides the correct spectral 
index. We found that the disagreement is mainly due to the steeper tail and the overall normalization, which is overestimated by the 
model. A simple power-law with an exp onential cutoff, oc (Eq / E) a+1 (1 — e~ kT ' B ), where a is determined by the transcendental 
equation (Titarchuk & Lyubarskij 1995), gives a reasonable agreement with simulations up to ~ 3 MeV. The chosen normalization 
provides t 



le correct value of the amplification factor eq. ( B4 ) 



A. 3. Annihilation 



The emissivity of a thermal plasma due to the electron-positron annihilation is (Dermer 1984) 



S a (e,kT) 



kTKl(l/Q) 



(2x z +l) 



chv(7r-l)e 2 - e cr a (7 r ), 



(A3) 



where K n is the modified Bessel function of the second kind and of order n, 7 r i s the relative Lorentz factor of the colliding particles 
(invariant), and er a (7r) is the annihilation cross section ( Jauch & Rohrlich 1976 ). 

The near-Earth intensity of the narrow annihilation line from the disk plane can be estimated by the assumption that all positrons 
which hit the disk annihilate in it (two annihilation photons per positron) 



la 



n+c R\ 
4 D 2 



(A4) 



where n + is the number density of positrons in the outer corona and \n + c is the flux density t oward the disk surfac e, Rd is the disk 
radius, D = 2.5 kpc is the distance, and id is the inclination angle of the disk plane (id ~ 40°, Liang & Nolan 1984). 



B. COOLING OF ELECTRONS 

The electron cooling in a thermal plasma at low number density and small optical depth is not very effective. The main channels 
are: bremsstrahlung, Comptonization, and Coulomb interactions with ions (mainly protons). 

B.l. Bremsstrahlung 

For a pure hydrogen plasma the ep-bremssrahlung luminosity dominates the ee-bremsstrahlung luminosity in the non-relativistic 
limit while at relativistic energies, O p> 0.5, the ee-bremsstrahlung dominates. The total energy emitted per unit volume of plasma 
electrons by ep- plus e + e~ -bremsstrahlung in the non-relativistic limit, <C 1, is (Haug 1985) 



,.NR , „,,NR ^ 2 3 , HZ I 1 



< P + <+ e - » ^a f rim e c*n;Ve + 2Z) +Z(1 + Z)\, (Bl) 
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where a; — 1/137 is the fine structure constant, and r e is the classical el ectron radius. The total ene rgy emitted by ee- plus 
e + e~-bremsstrahlung at the extreme-relativistic energies, 9 J> 1, is given by ( |Alexanian 1968j |Haug 1985| ) 



,.ER 



2Aa f rim e c 6 n z p (l + 2Z) 2 e{ln(26) - 0.5772 + 5/4}, 



and that for ep-bremsstrahlung is (von Stickforth 1961, paug 1975 ) 

= 12a f r 2 e m e c 3 n 2 p (l + 2Z)9{ln(29) - 0.5772 + 3/2}. 



(B2) 



(B3) 



B.2. Compton cooling 

An expression for the total energy losses of a plasma volume via Comptonization has been given by Dermer, Liang, & Canfield 
(1991) 



W Cr . 



L 



soft " 



1 - PA 



1 



-l-ln P/ InA 



(B4) 



where 



P 
A 



1 + 49 



K 3 (l/Q) 

K 2 (i/ey 



(B5) 



xq = Eo/kT (see eq. [ A2|), and L BO f t is the luminosity of the soft photon source. For r«l and io < 1 eq. (B4) is almost exact. 
The luminosity enhancement factor is given by r\ = L/L so f t = 1 + Wcm/ L so ft- 



B.3. Coulomb Coupling with Protons 

Stepney & Guilbert (1983) derived a general expression for the rate of energy transfer between populations of protons and electrons 
with Maxwellian distributions 



w C i 



\cn 2 p {l + 2Z) InA 



kTp k'J- 'fi 



K 2 (i/e e )K 2 (i/e p ) 



2K 



9 e + 9 p \ 2(9 e + 9 p ) 2 + l /9 e + 9 ? 



9 e 9 p 



9 e + 9 P 



9 e 9 p 



(B6) 



where InA is the Coulomb logarithm, 9 e = kT e /m e c 2 , and 9 P = kT p /m p c 2 are the dimensionless electron and proton temper- 
atures. The expression is symmetrical with respect to the electron and proton temperatures. In the limit of cold protons, 9 P — > 0, 
eq. (B6) reduces to 



w C i = ^—-Kric 6 ni(l + 2Z) InA 



26j + 29 e + 1 

K 2 (i/e e ) ' 



(B7) 



C. TRANSPORT CROSS SECTION 

Scattering at very small angles dominates in the Coulomb cross section, which reflects the long-range nature of the Coulomb 
interaction. However, for the diffusion process in plasma, small scattering angles are not very important. Additionally the ep- 
collisions are of minor importance compared to ee-collisions. We therefore restrict ourselves by considering the transport cross 
section only for ee-scattering, which provides us with an estimate on typical values of the relevant cross sections. 

The transport cross section for a test electron is defined by 



vtrM = I d 3 P2 V 7 ' 1 /(pa) / dQ (1 - cos 9)% 
J 7172 J dU 

= I d 3 P2 ^EIf{p 2 ) [dSl* (l -co S 6)^, (CI) 
J 7172 J ail* 

where f3i, 71 are the dimensionless speed and the Lorentz factor of the test particle, f(p 2 ) is the Maxwell-Boltzmann distribution, 
p 2 = 72/^2 is the momentum of the plasma particles, da/dfl is the differential cross section of the Coulomb scattering (J auch & 
Rohrlich 1976), and the asterisk marks the center-of-mass system (CMS) variables. The scattering angle, expressed in CMS variables, 
is cos 8 — (f3 c + cos 9*) / (1 + (3 C cos 9*), where /3 C , 7 C = (71 + 72 ) / \J2 (j r + 1 ) are the speed and Lorentz factor of the CMS relative 
to the laboratory system. Changing to the integration variables 72 and j r we obtain 

^^ = 2^0^(1/0) 1 **-tf-l) 1/2 J d l2 a tr (^ c )e-^, (C2) 
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where 7 ± = 717,- (1 ± PiPr), and 



Gtr{lr,lc) 



d(cos6*) 



(3 C + cos 9* 



1 + /3 C cos 9* 



da* 



d{cos9*) 



(C3) 



is the transport cross section for scattering angles 9* greater than the limiting angle 9* L — ► 0, where 9* M = tt/2 for M0ller scattering 
and 0^ = 7r for Bhabha scattering. 
For M0ller scattering the expression is 



CT*r(7r,7c) 



( 7r -l) 2 (7r + l) 



(47c -1) In 



l + P c \ , 2(1 -&) 
1 + /3 C 



(In A + In V2) + 1 



7r + l 



5 ) ln(l + /? c ) - — +41n2 + l 



(C4) 



The mean free path of a test electron in a thermal plasma would be 



A( 7 i;G) 



0i 



n_o- tr (7i; 9) 



(C5) 



REFERENCES 

Alexanian, M. 1968, Phys. Rev., 165, 253 
Abramowicz, M. A., et al. 1995, ApJ, 438, L37 

Balucihska-Church, M., Belloni, T., Church, M. J., & Hasinger, G. 1995, A&A, 302, L5 

Barr, P., White, N., & Page, C. G. 1985, MNRAS, 216, 65p 

Beall, J. H., Knight, F. K., Smith, H. A., & Wood, K. S. 1984, ApJ, 284, 745 

Bisnovatyi-Kogan, G. S., & Lovelace, R. V. E. 1997, ApJ, 486, L43 

Bowyer, S., Byram, E. T., Chubb, T. A., & Friedman, M. 1965, Science, 147, 394 

Chen, X., et al. 1995, ApJ, 443, L61 

Cui, W., Zhang, S. N., Focke, W., & Swank, J. H. 1997, ApJ, 484, 383 
Dermer, C. D., Liang, E. P., & Canfield, E. 1991, ApJ, 369, 410 
Dermer, C. D. 1984, ApJ, 280, 328 
Dolan, J. F. 1992, ApJ, 384, 249 

Done, C, Mulchaey, J. S., Mushotzky, R. F, & Arnaud, K. A. 1992, ApJ, 395, 275 
Dove, J. B., Wilms, J., Maisack, M., & Begelman, M. C. 1997, ApJ, 487, 759 

Ebisawa, K., et al. 1992, Frontiers of X-Ray Astronomy, Y. Tanaka and K. Koyama (Tokyo: University Academy Press), 301 
Ebisawa, K., et al. 1996, ApJ, 467, 419 

Esin, A. A., McClintock, J. E., & Narayan, R. 1997, ApJ, 489, 867 

Gierlihski, M., et al. 1997, MNRAS, 288, 958 

Haardt, F, Done, C, Matt, G., & Fabian, A. C. 1993, ApJ, 411, L95 

Haug, E. 1975, Z. Naturforsch., 30a, 1546 

Haug, E. 1985, A&A, 148, 386 

Haug, E. 1987, A&A, 178, 292 

Hua, X.-M., & Titarchuk, L. 1995, ApJ, 449, 188 

Inoue, H. 1989, 23rd ESLAB Symp. on Two Topics in X-Ray Astronomy (Noordwijk: ESA), 2, 783 
Jauch, J., & Rohrlich, F. 1976, The Theory of Photons and Electrons (New York: Springer) 
Jourdain, E., & Roques, J. P. 1994, ApJ, 426, LI 1 
Kemp, J. C, et al. 1983, ApJ, 271, L65 
Kitamoto, S., et al. 1990, PASJ, 42, 85 

Kolykhalov, P. I., & Sunyaev, R. A. 1979, Sov. Astron., 23, 189 

Kuznetsov S., et al. 1997, MNRAS, 292, 651 

Li, H., Kusunose, M., & Liang, E. P. 1996, ApJ, 460, L29 

Li, H., & Miller, J. A. 1997, ApJ, 478, L67 

Liang, E. P., & Dermer, C. D. 1988, ApJ, 325, L39 

Liang, E. P., & Nolan, P. L. 1984, Spa. Sci. Rev., 38, 353 

Liang, E. P. 1990, A&A, 227, 447 

Ling, J. C, Mahoney, W. A., Wheaton, Wm. A., & Jacobson, A. S. 1987, ApJ, 321, LI 17 
Ling, J. C, & Wheaton, Wm. A. 1989, ApJ, 343, L57 
Ling, J. C, et al. 1997, ApJS, 484, 375 

McConnell, M., et al. 1997, in AIP Conf. Proc. 410, Fourth Compton Symposium, ed. C. D. Dermer, M. S. Strickman, & J. D. 
Kurfess (New York: AIP), 829 



11 



MOSKALENKO, COLLMAR, & SCHONFELDER 



Vol. 999 



Moskalenko, I. V., Collmar, W., & Schonfelder, V. 1997, in AIP Conf. Proc. 4 10, Fourth Compton Symposium, ed. C. D. Dermer, 

M. S. Strickman, & J. D. Kurfess (New York: AIP), 863 ( |astro-ph/9709179| ) 
Moskalenko, I. V., & Jourdain, E. 1997a, in AIP Conf. Proc. 4 10, Fourth Compton Symposium, ed. C. D. Dermer, M. S. Strickman, 

& J. D. Kurfess (New York: AIP), 88 1 ( |astro-ph/9706283l) 

Moskalenko, I. V., & Jourdain, E. 1997b, A&A, 325,401 ( hstro-ph/970207l| ) 
Narayan, R., & Yi, I. 1995, ApJ, 452, 710 
Oda, M. 1977, Spa. Sci. Rev., 20, 757 

Owens, A., & McConnell, M. L. 1992, Comments Astrophys., 16, 205 

Phlips, B. F, et al. 1996, ApJ, 465, 907 

Priedhorsky, W., Terrell, J., & Holt, S. S. 1983, ApJ, 270, 233 

Pringle, J. E. 1981, ARA&A, 19, 137 

Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337 

Skibo, J. G., Dermer, C. D., Ramaty, R., & McKinley, J. M. 1995, ApJ, 446, 86 

Skibo, J. G., & Dermer, C. D. 1995, ApJ, 455, L25 

Stepney, S., & Guilbert, P. W. 1983, MNRAS, 204, 1269 

von Stickforth, J. 1961, Z. Physik, 164, 1 

Sunyaev, R. A., & Titarchuk, L. G. 1980, A&A, 86, 121 

Svensson, R. 1982, ApJ, 258, 335 

Svensson, R. 1984, MNRAS, 209, 175 

Tanaka, Y. 1991, Lect. Notes in Phys. 385, Iron Line Diagnostics in X-Ray Sources, A. Treves (Berlin: Springer), 98 

Titarchuk, L. G. 1994, ApJ, 434, 570 

Titarchuk, L., & Lyubarskij, Yu. 1995, ApJ, 450, 876 

Titarchuk, L. G, & Zannias, T, 1998, ApJ, 493, 863 



